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Abstract 

The Berezinskii-Kosterhtz-Thouless (BKT) mechanism describes the breakdown of superfluidity 
in a two-dimensional Bose gas or a two-dimensional gas of paired fermions. In the latter case, a 
population imbalance between the two pairing partners in the Fermi mixture is known to influence 
pairing characteristics. Here, we investigate the effects of imbalance on the two-dimensional BKT 
superfluid transition, and show that superfluidity is even more sensitive to imbalance than for three 
dimensional systems. Finite-temperature phase diagrams are derived using the functional integral 
formalism in combination with a hydrodynamic action functional for the phase fluctuations. This 
allow to identify a phase separation region and tricritical points due to imbalance. In contrast 
to superfluidity in the three-dimensional case, the effect of imbalance is also pronounced in the 
strong-coupling regime. 
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I. INTRODUCTION 



Quantum phenomena which occur at very low temperatures are a subject of intense 
experimental and theoretical study. Recent progress in the experimental investigation of 
ultra-cold atoms stimulated an unprecedented interest to the theoretical problems of con- 
densation of cold bosons and pairing of interacting fermions (see, e. g., the review |l| and 
references therein). These phenomena are related to a variety of objects including stars, 
dense nuclear and quark matter, and plasma systems 

Phase transitions of quantum systems strongly depend on their dimensionality. Two- 
dimensional Fermi gases have remarkable features, which are not observed in three di- 
mensions. The quasi-2D regime for cold atoms can be reached using a sufficiently strong 
confinement of atoms along one direction, so that they occupy only the lowest size quan- 
tization subband. Advances in pairing of trapped cold atomic Fermi gases with a con- 
trolled geometry of a trapping potential allow experimentalists to realize systems of different 
(quasi) dimensionalities. Both condensation of bosonic atoms [sl and pairing of fermions 
has been observed recently in optical lattices. 

Q 

The Mermin-Wagner-Hohenberg theorem [7| shows that in a uniform, two-dimensional 
(2D) system, the long range order is destroyed by thermal fluctuations so that 2D Bose gases 
cannot undergo Bose-Einstein condensation at nonzero temperatures 0,[9|,[l^. Neverthe^ 
these two-dimensional systems can form a " quasicondensate" and exhibit superfluidity 
Kosterlitz and Thouless 12] showed that the mechanism driving the superfluid-to-normal 
phase transition in this case is the proliferation of vortices and antivortices above a critical 
temperature Tbkt, spoiling the phase coherence. The experimental situation differs in 
that considered by Mermin, Wagner and Hohenberg in that the gas is not uniform. The 
presence of a trapping potential in the experimental realizations changes the density of states, 



ess, 



111. 



and in a conflned 2D system, true Bose-Einstein condensation does become possible 



13|. 



However, the Berezinskii-Kosterlitz- Thouless (BKT) mechanism continues to play a role in 

n 

the suppression of superfluidity also in trapped systems|14l]. Indeed, a BKT-type crossover 
in a trapped quantum degenerate gas of ^''Rb has recently been observed by Hadzibabic 



et al. 



15| using an interference technique to detect phase defects |16l|. This has stimulated 



renewed theoretical interest in the BKT transition in trapped dilute gasesjlTj. 

The discussion in the above paragraph focuses on Bose gases. Rice [l8|, and Schmitt- 
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Rink, Varma, and Ruckenstein [19| extended these results to the case of a gas of paired 
fermions. The pairing characteristics depend on the the interatomic interaction strength. 
The su per fluid phase transition in a 2D Fermi gas at T = has been studied by Randeria 
et al. 2m[ in a wide range of interaction strengths from the weak-couphng hmit, where 
the Bardeen-Cooper-Schriffer (BCS) pairing regime is reahzed, to the strong-couphng hmit 
which corresponds to the Rose condensation of interacting bound pairs. Petrov et ah 21 1 
analytically investigated pairing of a quasi 2D Fermi gas at a finite temperature in the weak- 
coupling BCS limit and in the strong-coupling limit where the BKT superfluid transition 
occurs. 

The influence of spin imbalance on the superfluid properties of a Fermi gas is of particular 
interest. A phase separation between the superfluid and normal component of an interact ing 



Fermi gas of cold atoms with unequal spin populations in 3D was observed in Refs. 



22 



23|, 



24l |. In these experiments, both the interaction strength between ultra-cold fermions and the 



population imbalance can be controlled using the Feshbach resonance. Allowing the control 
of population or mass imbalance opens a unique possibility to investigate the stability of 
fermion pairing and, in particular, to determine experimentally the equation of state for an 
imbalanced fermion system |24|. 

The BKT phase transition for a balanced 2D Fermi gas with the s-wave scattering in 
the balanced case has been theoretically analyzed by Botelho and de Melo {25 1. They 
treated the fermion pairing in 2D by the path integral technique taking into account phase 
fluctuations. Within a similar approach, using an effective Hamiltonian which involves 
fermions interacting with each other and with dressed molecules, the BKT transition has 
been considered for a quasi-2D trapped Fermi gas |26 |. 



In the present work, we extend the approach of Ref. [25j to investigate the effect of 
population imbalance on the BKT phase transition in a 2D Fermi gas. Using the Hubbard- 
Stratonovich transformation, we derive a hydrodynamic effective bosonic action, which in 
the limiting case of a balanced gas is reduced to the effective action of Ref. [251]. On the basis 
of the obtained effective action, we analyze phase diagrams for an interacting, imbalanced 
Fermi gas in 2D. The paper is organized as follows. In Sec. [Tll we describe the theoretical 
formalism for interacting fermions in 2D. In Sec. IIIIl we analyze the dependence of the 
critical temperature of the BKT transition on the coupling strength and on the population 
imbalance. The section is followed by conclusions. Sec. HVl 
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II. FUNCTIONAL INTEGRAL DESCRIPTION 



A. General formalism 



We consider a two- component gas of interacting fermions in 2D, with the s-wave pairing 
and with a population imbalance. The partition function of the system of fermions in 2D is 
expressed as the path integral over Grassmann variables ['?/'x,t,(t, '^x,t,(t] , 



Z (X j V [V^x,<7 (r) , V^x,,T (r)] exp (-5) . 



The action functional 



S — Sq + Sint 

is a sum of the free-fermion and interaction contributions, 



Sq = J dr j d^x 



S; 



int 



2„ / j2 



^x,<T (^) ( ^ - Vx - /ia ) i^^,a (r) 



dr d X d yV{x-y) ^x,t (^) ^y,i (^) ^y,i (^) V-x,! (^) 



(1) 



(2) 



(3) 



(4) 



where f3 = is the inverse to the temperature. We use the units in which h = 1, the 
fermion mass m = 1/2, and the Fermi energy Ef = {^-nnof^ I (2m) = 1 (where no is the 
fermion density in 2D). We express the results below in terms of the averaged chemical 
potential /i = (/i| + ^i) /2 determining the total number of fermions and the imbalance 
potential Q = {^^ — /i i ) / 2. For the interaction potential, we use the separable expression 



proposed in Refs. [25|, l27| , 



^k,k' = ^rkFk', (5) 
where g is the interaction strength. The factor Fk describes a finite-range potential, 

k N 



1 + 



kn 



(6) 



where R ~ k^-^ plays the role of the interaction range. The particular case of the contact 
interaction corresponds to /eq ^ oo so that Fk — 1. The interaction term Sint of the action 
functional is then given by 



S., 



int 



t) B„ (t) . 



(7) 
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where L is the lateral size of the 2D system, and the collective coordinates i?q (r) are 
determined as 

(^) = 5Z Tka-k+f ,i {r) ak+a,| (r) . 

k 

Further on, we apply the Hubbard-Stratonovich (HS) transformation. Introducing the ex- 
tended action 

Sext = S ~ 



r dT-J2kir)^^{T) (8) 
-^0 9 „ 



with the auxiliary Bose field (HS field) 0q (r), and performing the shift of boson coordinates, 
which eliminates the fermion-fermion interaction term Sint, the HS transformation results 
in the action 

r-/3 



Sext = So- [ dTjYl (^) (^) + (^) (^)] 

Jo ^ q 

dr- Yl (^) '^'i (^) • (9) 



■/3 

' ^ 



q 



Because the phase fluctuations about the saddle point are, in general, not small, the 



boson (HS) and fermion variables in the coordinate representation are transformed as 28| : 



0r (r) = 0; (r) e^^^(-), V^r,. (r) = V^;,. (r) e^^-(-)/2. (10) 

In the same formalism for a Fermi gas in 3D, the further step is the path integration over 
fermion variables and the expansion of the resulting bosonic action over fluctuations about 



the saddle point [291]. This method provides a description of a superfluid phase transition 
between the normal phase and the true condensate of fermion pairs. In the 2D case, at least 
when restricting the expansion by quadratic fluctuations, the superfluidity occurs only at 



T = 19| . Traven [30] showed that the interaction between fluctuations of the pairing field 
in a 2D attractive Fermi gas allows a superfluid phase transition at a very low temperature. 
However, the superfluid state can exist in a 2D Fermi gas at relatively high temperatures 
as a quasicondensate, i. e., a state with fermion pairing "where the density fluctuations are 



suppressed but the phase still fluctuates" [3l|. In two dimensions, the quasicondensate can 
be realized through bound vortex-antivortex pairs 
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B. Fluctuating phase 



After integrating out the fermion (Grassmann) variables, an affective action in the bose 
field 0r is obtained. The remaining functional integral over this Bose field cannot be taken 
in general. Several levels of approximation can be made to get results. The crudest approx- 
imation is the mean-field approximation which replaces the field by a constant, 0r \^\, 
the saddle point. To improve on this, fluctuations around the saddle point can still be taken 
into account; this can be done in an exact way only up to quadratic order in the fluctu- 
ation. One can choose to write the fluctuations as 0r —* |A| + 6r, with 6r complex; this 
corresponds to the NSR approach that investigates the presence of a real condensate rather 
than a quasicondensate. Alternatively the NSR fluctuations can be written as amplitude 



and phase fluctuations. 



A I + \5r\) e where \6r\ and e " are real fields. In order to 



treat the quasicondensate state, relevant for the 2D case, we have to focus on phase fluctu- 



ations: that corresponds to setting 



lA 



I A| e*^"'. We will simplify the notation and write 



A for the energy gap parameter of the Bogoliubov excitations. Moreover, we will 
assume that the remaining fluctuation field 6r varies slowly as a function of position and 
time with respect to the variations of the fermion fields. A similar assumption was used in 
28 1 . The resulting hydrodynamic action is structurally similar to the saddle-point 



Refs. 



action for imbalanced fermions in 2D 



32| 



S, 



eff 



- / dr 







- In ( 2 cosh /5Ck + 2 cosh /3Ek 

P 



in which the fermion energy = /c^ — /x, the Bogolubov excitation energy = \Ji'^ + A^F^ 
and the imbalance potential ^ are replaced by expressions depending on boson coordinates: 



Ck = c - v^^ • k 



i do 1 ,„„,2 i 



(12) 
(13) 
(14) 



Keeping the quadratic-order fluctuation terms we arrive at the action Se/y as the sum 
of the saddle-point action S'^p, which coincides with the expression (5) of Ref. |32|, and the 
diagonal quadratic form of the phase fluctuations 



S 



fi 



/3 



dr d r 



(15) 
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The coefficients in the fluctuation action are the pair superfluid density 

in/ c. \ 1.2 



2 V y 2 



(16) 



and the constant 



where Xu the function 



A 



—y 



El 



Xk + -^Yk 



El 



sinh 13 Ek 



cosh /3C + cosh (3Ek 



and Yk is the extension of the Yoshida distribution to imbalanced fermions: 



Yk = f3 



cosh pC, cosh PEk + I 



(17) 



(cosh /3( + cosh (3Ek) 

The action in Eq. fllSp describes a Bose gas of spin waves [25f] with the energy spectrum 
ujk given by 

Uk = ck, c= -y/^. (18) 
The spin wave contribution to the thermodynamic potential is given by the expression 

(19) 



l5:in(l-e-/^-) 



The total thermodynamic potential of the fermion system taking into account phase fluc- 
tuations is the sum of the spin-wave term (JT9l) and the saddle-point contribution, which is 
provided by the saddle-point action, 



sp 



/3 



In (2 cosh /3C + 2 cosh pE^) - 
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(20) 



For an imbalanced quasi 2D Fermi gas in an optical potential, the mean-field zero- 
temperature phase diagrams were analyzed in Ref. 32| on the basis of this saddle-point 
action (neglecting spin- wave contributions). 

The gap parameter A for an imbalanced Fermi gas is determined through the minimiza- 
tion of the saddle-point thermodynamic potential Qsp as a function of A at given /?, /i, (: 



'^sp 



dA 



0. 



(21) 



For a complete determination of thermodynamic parameters at a given temperature, the 
minimum condition (1211) is solved jointly with the number equations: 

where n and 6n are the total fermion density n and the density difference 6n, respectively. 



C. BKT transition temperature 



The coupled gap and number equations (^T\i and (1^ have to be solved together. In 
equation ( l22l) . three different levels of approximation can be made, in analogy to the ap- 
proximations for the bose field 0r(T) as discussed in the beginning of the previous subsection. 



The first and simplest case is the mean-field approximation (as described in Refs. |25l. |26| ) . 
where we use fl = Qgp- This corresponds to setting the (Hubbard-Stratonovic) Bose field 
equal to a constant (p = Amf (both constant in amplitude and in phase). The constant can 
be determined by minimizing Qgp and allows to determine the phase transition line between 
the normal state, in which Amf = 0, and the quasicondensate of fermion pairs, in which 
Amf 7^ 0. The temperature separating the aforementioned phases will be denoted by Tmf- 
For a balanced gas, the phase transition between normal and paired states is of the second 
order. For an imbalanced gas, also the first-order phase transition between the normal and 
paired states is possible. 

The paired state below Tmf is not always superfluid. Proliferating vortices and antivor- 
tices destroy phase coherence in the quasicondensate, and suppress superfluidity [l^ . The 
mean-field approximation does not allow to determine the temperature Tbkt separating the 
superfluid quasicondensate state from the non-superfiuid paired state. To investigate the 
question of superfiuidity and determine the BKT temperature, we need to use Q = Qgp + ^sw 
in 0221) . where Qgw is the free energy contribution 0191) from the phase fiuctuations. This 
corresponds to giving the (Hubbard-Stratonovic) Bose field a constant amplitude, but al- 



lowing its phase to fiuctuate freely, 
universal Nelson-Kosterlitz equation 



he BKT transition temperature is then a root of the 



33| 



Tbkt - ^Ps {Tbkt) = 0. (23) 

As distinct from the aforesaid phase transition at T = Tmf, the BKT phase transition at 
T = Tbkt is characterized by an abrupt change of the superfiuid density from zero to a 
finite value satisfying Eq. 016p . The phase- fiuctuation contribution to the density vanishes 
at the phase boundary between the paired state and the normal state, because at Amf = 
the superfiuid density ps turns to zero. 

Finally, the approximation of setting Q = Qgp + ^sw in ( l22l) does not allow to determine 
whether the superfluid quasicondensate can form a true condensate below some critical 
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temperature Tbec- Note that the emergence of a true condensate is only possible in a finite 
system, as Tbsc ^ for L ^ oo. To determine the critical temperature we should use 
VL = Qsp + ^fiucts in (l22l) . where ^l/iucts contains contributions from both amplitude and 
phase fiuctuations in the Hubbard- Stratonovic Bose field. 

III. RESULTS AND DISCUSSION 

When Tmf > 0, equation (123!) can be satisfied if the superfiuid density ps 7^ 0. At 
T = Tmf, the superfiuid density for a balanced Fermi gas turns to zero. Therefore for a 
balanced Fermi gas Tbkt is always lower than Tmf- On the contrary, in the case of unequal 
"spin up" and "spin down" populations, for a sufficiently high population imbalance the 
phase transition at T = Tmf can be of the first order, when the gap parameter A changes 
discontinuously from zero to a finite value. In this shown below, Tbkt and Tmf 

can coincide in a definite range of the coupling strength. 

In the region between Tbkt and Tmfi the superfiuidity is destroyed owing to a prolifer- 
ation of free vortices. However, the phase of a 2D Fermi gas in this region is not a uniform 
normal state, because the gap parameter for Tbkt < T < Tmf can be other than zero, so 
that fermion pairing can occur in that temperature region. The region between Tbkt and Tc 
can be therefore attributed to a state in which pairing can occur but the phase coherence is 
destroyed. For slightly higher temperatures with respect to Tbkt, vortices form a disordered 
gas of phase defects jl| . For higher temperatures the concept of vortices is inapplicable due 
to density fiuctuations. 

In general, there are the following regions in phase diagrams for an imbalanced Fermi 
gas in 2D: (1) the normal state in which the gap parameter A = 0, (2) the superfiuid 
quasicondensate state in which A 7^ and ps 7^ 0, (3) the state in which pairing is possible 
but without phase coherence, and (4) the phase separation region, which appears due to 
a population imbalance. In the latter one, no solution exist for the set of equations (121 p 
and ([22]) • Therefore, in the phase separation region the system will unmix in a phase with 
lower (or no) imbalance and the normal phase at higher imbalance. At the first-order phase 
transition, the system abruptly passes through the phase separation region. 

In Fig. 1, we represent phase diagrams in the variables {T, Eh}, where T is the tempera- 
ture, and Ef, is the binding energy of a two-particle bound state in 2D. The phase diagrams 
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in Fig. 1 are obtained at a given value of the imbalance chemical potential (. The top 
panel shows the balanced case (C = 0), and the middle and lower panels illustrate how the 
diagram changes with increasing difference between the chemical potentials of up and down 
species. The energy Eb is the parameter which characterizes the coupling strength of the 
fermion-fermion attractive interaction. In two dimensions, the strength g of the contact 



interaction is related to this binding energy through 20[ |. 



1 m 



In {E/Ei 



b 



(24) 



The binding energy itself can be related to the experimental parameters. The two- 
dimensional system is created through a strong confinement of the third direction; in general 
this strong confinement can be associated with a harmonic potential with frequency uj^ (and 
oscillator length II). The two-particle bound state exists for all values of the (3D) s-wave 
scattering length of the fermionic atoms and its binding energy is given by: 



Chhjj 

TT 



i^, = r:r^£exp v^^), (25) 



with C ^ 0.915 (cf. Ref. [SJ). 

The dashed curves in Fig. 1 correspond to the phase transition curves, TMpiEt), within 
the mean-field approach. When the system is imbalanced, two changes occur with respect 
to the balanced case: (i) below a certain value of the binding energy, Tmf 0, and (ii) a 
tricritical point appears on this curve. The mean-field phase transition curve splits below 
this tricritical point and the phase separation region opens up. In this tricritical point, there 
is coexistence of the paired state without phase coherence, the normal state, and the phase 
separation region. The area of the phase separation region is narrow and broadens with 
increasing the imbalance potential. 

The full curve corresponds to the BKT phase transition line TsKriEb)- For the bal- 



anced case, we retrieve the results obtained in Ref. j25|] and find that the BKT transition 
temperature is nonzero for any finite binding energy. In particular, we find that in the 
strong-coupling limit (for large Ef,), the superfluid density ps rises to its maximum, equal to 
half the total fermion density, so that the BKT transition temperature obtained from Eq. 
( I23l) tends to a finite value T^^^^ = (l/8)Tp. When imbalance is introduced, two qualitative 
changes occur in the BKT transition curve. Firstly, below a critical value of the binding 
energy, Tbkt 0. Imbalance suppresses BKT-superfluidity as it does for 3D-superfluidity. 
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The critical binding energy needed to restore superfluidity shifts to higher values as the im- 
balance potential grows. Secondly, there appears a tricritical point also on the BKT curve, 
indicating coexistence of the superfluid state, the state with pairing but no phase coherence, 
and the phase separation region. 

The dotted curve in the region of the normal state separates two regimes: (i) the regime 
where the thermodynamic potential contains only one minimum at A = above the dotted 
curve, and (ii) the regime where the thermodynamic potential contains two minima: the 
lower minimum at A = and a higher one at a value A 7^ 0. This second minimum suggests 
the existence of a metastable superfluid state. 

In Fig. 2, we choose to flx 5n/n, the relative population imbalance itself, rather than the 
imbalance potential and study the phase diagram in the variables {Ef,, T, The inset 
shows a slice of this phase diagram at Ef, = 0.01. At this binding energy, the BKT transition 
temperature is a non-single- valued function of 6n/n. This result is not unexpected, because 
a similar non-single- valued behavior of the critical temperature exists also for the superfluid 
phase transition in a 3D Fermi gas jssl. In the 3D case this can be related to the Sarma 



mechanism [36|]: at nonzero temperatures, a balanced superfluid system coexists with an 
imbalanced gas of thermal excitations - these excitations carry some of the excess majority 
component of the imbalanced gas. 

The BKT phase transition for a Fermi gas in 2D is suppressed at lower values of the pop- 
ulation imbalance than the superfluid phase transition for a Fermi gas in 3D. Furthermore, 
this suppression is more strongly revealed at higher binding energies. The reason of such 
a dependence of the BKT phase transition as a function of the binding energy consists in 
the following. Let us consider flrst a part of the phase diagram in the variables {Eb,T,^) 
taking into attention only the mean-fleld phase transition. Similar phase diagrams were 
calculated for a 3D Fermi gas in Ref. 35] . An important feature of phase diagrams for both 
3D and 2D imbalanced Fermi gases is the fact that the phase separation region lies at lower 
temperatures with respect to the superfluid state. This means that, when the temperature is 
gradually lowered at a flxed imbalance 6n/n (which is assumed to be sufficiently low to pass 
the region of the superfluid state), we can observe the superfluid phase transition from the 
normal state to the superfluid state, and then, when the temperature is lowered further, the 
superfluid state becomes non-stable in the phase separation region. So, a sufficient lowering 
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of the temperature leads to an instability of the superfluid state of an imbalanced Fermi 
gas. On the one hand, with increasing the coupling strength, the phase separation region 
broadens, and its upper-temperature boundary shifts to higher temperatures. On the other 
hand, the BKT transition temperature does not unrestrictedly increase with an increasing 
coupling strength. As follows from the Nelson-Kosterlitz equation ( |23l) . the upper bound 
for Tbkt is Tbkt = (l/8)Tp. This is the case when the superfluid pair density achieves its 
maximal possible value for a balanced gas pi™^^"* = = As a result, with an increasing 
coupling strength, we arrive at a stronger overlap between the area of the superfluid quasi- 
condensate state (below Tbkt) and the phase separation region. It is natural to expect that 
at sufficiently high binding energy, the phase separation region can completely cover the area 
below Tbkt, so that the superfluid state vanishes at all temperatures. In this connection, 
the study of phase diagrams for a 2D Fermi gas in the variables (T, Eh) at a given relative 
population imbalance (i. e., slices of the above considered three-dimensional phase diagram 
at fixed ^) represents a particular interest. 

In Fig. 3, we have plotted phase diagrams for a 2D Fermi gas in the variables (T, Eb) at 
the relative population imbalance ^ = 0.03 and 0.1. We see that even at a low population 
imbalance — = 0.03, the region where the superfiuid state exists is relatively narrow. Fur- 
thermore, with an increasing binding energy E^, the upper-temperature bound of the phase 
separation region rises faster than the BKT transition temperature. As a result, the range of 
coupling strengths at which the BKT superfiuid phase transition can occur, is restricted both 
from lower and higher coupling. For the comparison, the regular superfiuid phase transition 
in a 3D imbalanced Fermi gas is restricted only from lower couplings js^]. We can see from 
Fig. 3 (6) that for a higher (but relatively low) imbalance ^ = 0.1, the higher-temperature 
bound of the phase separation region lies higher than Tbkt for all considered values of the 
binding energy E^. As a result, for ^ = 0.1 the superfiuid state is absent. The obtained 
behavior of phase diagrams for a 2D Fermi gas confirms our suggestions above. 

IV. CONCLUSIONS 

In summary, we have described the effects of imbalance on the Berezinskii-Kosterlitz- 
Thouless superfluid transition in a 2D Fermi gas, through the functional integral formalism. 
Owing to a population imbalance, the superfluid state cannot exist for all values of the 
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coupling strength, but only above a certain critical binding energy which depends on the 
imbalance. The larger the imbalance potential, the higher is this critical binding energy. 
As distinct from the balanced case, there is a phase-separated region in phase diagrams, 
in which no uniform state can exist. As a result, tricritical points appear at the phase 
diagrams, in which three different regimes coexist. The area of the phase-separated region 
increases with increasing the population imbalance. Due to the rise of the upper-temperature 
bound of the instability region, the area of the superfiuid state at a fixed relative population 
imbalance decreases with increasing the binding energy. Therefore a population imbalance is 
a factor destroying superfluidity in 2D systems, especially at high binding energies. The BKT 
transition can be experimentally observable for an imbalanced Fermi gas through the phase 
separation in a quasi 2D trap. As follows from the obtained results, the parameters of the 
state of that system (e. g., critical temperatures and/or density profiles) are expected to be 
substantially more sensitive to the population imbalance than the corresponding parameters 
for a 3D gas. 

Acknowledgments 

The authors would like to acknowledge fruitful discussions with M. Wouters, M. 
Oberthaler, and H.T.C. Stoof. This work was supported by FWO-V projects G. 01 15.06, 
G.0356.06, G.0180.09N, G.0370.09N, the WOG WO.033.09N (Belgium), and INTAS Project 
no. 05-104-7656. J.T. gratefully acknowledges support of the Special Research Fund of the 
University of Antwerp, EOF NOT UA 2004. 



[1] 1. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008). 

[2] R. Casalbuoni and G. Nardulli, Rev. Mod. Phys. 76, 263 (2004); Yu.A. Litvinov et al, Phys. 

Rev. Lett. 95, 042501 (2005); E. Braaten and H.-W. Hammer, Phys. Rep. 428, 259 (2006); 

G. Labeyrie, F. Michaud, and R. Kaiser, Phys. Rev. Lett. 96, 023003 (2006). 
[3] P. F. Bedaque, Nud.Phys. A 697, 569 (2002); M. G. Afford, J. A. Bowers and K. Rajagopal, 

Phys.Rev. D 63, 074016 (2001). 
[4] A. Sedrakian and U. Lombardo, Phys. Rev. Lett. 84, 602 (2000). 

[5] M. Greiner, M. O. Mandel, T. Esslinger, T. Hansch, and I. Bloch, Nature 415, 39 (2002). 

13 



[6] J. K. Chin, D. E. Miller, Y. Liu, C. Stan, W. Setiawan, C. Sanner, K. Xu, and W. Ketterle, 

Nature (London) 443, 961 (2006). 
[7] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966); P. C. Hohenberg, Phys. Rev. 

158, 383 (1967). 

[8] N. N. Bogolubov, Lectures on Quantum Statistics (Gordon and Breach, New York, 1970), Vol. 
2. 

[9] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics (Pergamon Press, Oxford, 1980), Pt. 2. 
[10] V. N. Popov, Functional Integrals in Quantum Field Theory and Statistical Physics (D. Reidel 

Publishing, Dordrecht, 1983). 
[11] V. L. Berezinskii, Sov. Phys. JETP 32, 493 (1971). 

[12] J. M. Kosterlitz and D. J. Thouless, J. Phys. C 6, 1181 (1973); J. M. Kosterlitz, J. Phys. C 
7, 1046 (1974). 

[13] D. S. Petrov, M. Holzmann, and G. V. Shlyapnikov, Phys. Rev. Lett. 84, 2551 (2000); V. 
Bagnato, D. Kleppner, Phys. Rev. A 44, 7439-7441 (1991); J. Tempere and J. T. Devreese, 
Solid State Comm. 101, 657 (1997); J.P. Fernandez and W.J. Mullin, J. Low Temp. Phys. 
128, 233 (2002). 

[14] N. Prokof'ev and B. Svistunov, Phys. Rev. A 66, 043608 (2002); S. Sachdev and E. Demler, 

Phys. Rev. B 69, 144504 (2004). 
[15] Z. Hadzibabic, P. Kriiger, M. Cheneau, B. Battelier, and J. Dalibard, Nature 441, 1118 (2006); 

P. Kriiger, Z. Hadzibabic, and J. Dahbard, Phys. Rev. Lett. 99, 040402 (2007). 
[16] Z. Hadzibabic, S. Stock, B. Battelier, V. Bretin, and J. Dalibard, Phys. Rev. Lett. 93, 180403 

(2004); S. Stock, Z. Hadzibabic, B. Battelier, M. Cheneau, and J. Dalibard, Phys. Rev. Lett. 

95, 190403 (2005). 

[17] M. Holzmann, G. Baym, J.-P. Blaizot, and F. Lalo6, PNAS 104, 1476 (2007); Z. Hadzibabic, 
P. Kriiger, M. Cheneau, S.P. Rath, J. Dahbard, New Journal of Physics 10, 045006 (2008); 
T. P. Simula, M. J. Davis, and P. B. Blakie, Phys. Rev. A 77, 023618 (2008); M. Holtzmann, 
W. Krauth, Phys. Rev. Lett. 100, 190402 (2008); R.N. Bisset, M.J. Davis, T.P. Simula, P.B. 
Blakie, Arxiv preprint [arXiv:0804.0286] (20081. 

[18] T. M. Rice, Phys. Rev. 140, A1889 (1965). 

[19] S. Schmitt-Rink, C. M. Varma, and A. E. Ruckenstein, Phys. Rev. Lett. 63, 445 (1989). 
[20] M. Randeria, J. M. Duan, and L. Y. Shieh, Phys. Rev. B 41, 327 (1990). 

14 



[21] D. S. Petrov, M. A. Baranov, and G. V. Shlyapnikov, Phys. Rev. A 67, 031601 (2003). 

[22] M. W. Zwierlein, A. Schirotzek, C. H. Schunck, and W. Ketterle, Science 311, 492 (2006); 
M. W. Zwierlein, C. H. Schunck, A. Schirotzekl, and W. Ketterle, Nature (London) 442, 54 
(2006). 

[23] G. B. Partridge, W. Li, R. I. Kamar, Y.-A. Liao, and R. G. Hulet, Science 311, 503 (2006). 

[24] Y. Shin, C. H. Schunck, A. Schirotzek, and W. Ketterle, Nature (London) 451, 689 (2008). 

[25] S. S. Botelho and C. A. R. Sa dc Mclo, Phys. Rev. Lett. 96, 040404 (2006). 

[26] W. Zhang, G.-D. Lin, and L.-M. Duan, Phys. Rev. A 78, 043617 (2008). 

[27] R. D. Duncan and C. A. R. Sa de Melo, Phys. Rev. B 62, 9675 (2000). 

[28] S. De Palo, C. Castellani, C. Di Castro, and B. K. Chakraverty, Phys. Rev. B 60, 564 (1999). 

[29] C.A.R. Sa de Melo, M. Randeria, and J.R. Engelbrecht, Phys. Rev. Lett. 71, 3202 (1993). 

[30] S. V. Travcn, Phys. Rev. Lett. 73, 3451 (1994). 

[31] D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, Phys. Rev. Lett. 85, 3745 (2000). 

[32] J. Tempere, M. Wouters, and J. T. Devreese, Phys. Rev. B 75, 184526 (2007). 

[33] D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 1201 (1977). 

[34] D.S. Petrov and G.V. Shlyapnikov, Phys. Rev. A 64, 012706 (2001). 

[35] M. M. Parish, F. M. Marchetti, A. Lamacraft, B. D. Simons, Nature Physics 3, 124 (2007). 

[36] G. Sarma, J. Phys. Chem. Solids 24, 1029 (1963). 

[37] J. Tempere, S. N. Klimin, J. T. Devreese, and V. V. Moshchalkov, Phys. Rev. B 77, 134502 
(2008). 



15 



Figure captions 



Fig. 1. Phase diagrams for a 2D imbalanced Fermi gas in the variables (T, E^) for different 
values of the imbalance potential (. The white area is the normal state, the blue area is the 
superfluid state, the green area is the state with pairing and without phase coherence, and 
the grey area is the phase separation region. The dotted curve indicates the BKT transition 
to a metastable superfluid state. 

Fig. 2. Three-dimensional phase diagram for a 2D imbalanced Fermi gas in the variables 
(E^ft, T, ^) The slice at Ei, — 0.01 shows the phase diagram in the variables (T, ^) with the 
same denotations as in Fig. 1. 

Fig. 3. Phase diagram for an imbalanced 2D gas of cold fermions in the variables (T, Ei,) 
at the relative population imbalance ^ = 0.03 (a) and ^ = 0.1 (b). The denotations are 
the same as in Fig. 1. 



16 



Figures 



0.2 



a Mean-field / ^= 

/ Pairing witliout 



colierence 




Superfluid state 



0.0 
10 

0.3 



0.2 



0.1 



10- 



10 



10" 



/I ^ F 

Normal // ^ . . 

/ / Pairing 

/ / without 

/ / Phase coherence 

/ jq separation 

/ / 
/ / 





0.0 
0.2 

0.6 fc 



0.4 



0.2 




0.3 0.4 0.5 0.6 0.7 0.8 



Normal 
state 



Pairing 
without 
coherence 



Phase 



0.0^ 
0.6 



.-/■ i 

/separation/ Superfluid 



state 



0.8 1 1.2 1.4 1.6 1.8 2 

h F 



Fig. 1 J. Tempere et al. 



17 



E/E^= O.Ol —--Mean-field 
BKT 




Fig. 2 J. Tempere et al. 



18 



15 



E5; 



0.3 

0.2 

0.1 

0.0 
0.20 

0.15 

■^0.10 
0.05 



Mean-field 


/ 

/ a 


BKT 


5n/n = 0.03 


/ 


- 

Normal 


/ without 


state 




/' 




^ ■ 

^ ^^^^^^ 




Phase separation 



0.01 



0.1 



0.00^ 



Mean-field 


/'Pairing /' ^ 


BKT 


/ without / 


" hiln = 0. 1 , 


1 coherence^/^ 


/ 


/ 

✓ 


Normal / 


y 


state / 


y 

y 


/ 


Phases^aratim^/^^ 


0.01 


0.1 ] 




F 



Fig. 3 J. Tempere et al. 



19 



